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Abstract 

We develop an analytical theory for generic disorder-driven quantum phase transitions. We apply 
this formalism to the superconductor-insulator transition and we briefly discuss the applications to 
the order-disorder transition in quantum magnets. The effective spin-^ models for these transitions 
are solved in the cavity approximation which becomes exact on a Bethe lattice with large branching 
number iiT S> 1 and weak dimensionless coupling g <^ 1. The characteristic features of the low 
temperature phase is a large self-formed inhomogeneity of the order-parameter distribution near 
the critical point K > Kc{g) where the critical temperature Tc of the ordering transition vanishes. 
We find that the local probability distribution P{B) of the order parameter B has a long power-law 
tail in the region where B is much larger than its typical value Bq. Near the quantum critical point, 
aX K ^ Kc{g), the typical value of the order parameter vanishes exponentially, Bq oc e'^^^^"^"^^^^ 
while the spatial scale Ni^h of the order parameter inhomogeneities diverges as {K — Kc{g))~^. In 
the disordered regime, realized at K < Kc{g) we find actually two distinct phases characterized by 
different behavior of relaxation rates. The first phase exists in an intermediate range of K*{g) < 
K < Kc{g). It has two regimes of energies: at low excitation energies, uj < u!d{K, g), the many-body 
spectrum of the model is discrete, with zero level widths, while at u > ojd the level acquire a non- 
zero width which is self-generated by the many-body interactions. In this phase the spin model 
provides by itself an intrinsic thermal bath. Another phase is obtained at smaller K < K*{g), 
where all the eigenstates are discrete, corresponding to full many-body localization. These results 
provide an explanation for the activated behavior of the resistivity in amorphous materials on the 
insulating side near the SI transition and a semi-quantitative description of the scanning tunneling 
data on its superconductive side. 

PACS numbers: 03.67 
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I. INTRODUCTION 



Recently the subject of zero temperature quantum phase transitions and of the corre- 
sponding quantum critical points in translationally invariant systems got a lot of attention, 
the theoretical description of this phenomenon is mostly complete.- Much less is known 
and understood about the transition driven by the competition of a strong disorder and 
interactions in quantum systems which is the subject of this paper. 

The goal of the paper is twofold: to formulate a theoretical model that is relevant for 
the description of a number of experimental systems and to solve this model in the simplest 
controlled approximation. The physical systems that we shall focus on are disordered su- 
perconductors but the main results can also be applicable to disordered magnets, especially 
disordered ferromagnets in a random field. The new physics introduced by the strong disor- 
der is the appearance of new phases^ in which all or some excitations are localized in space 
and have infinite lifetime and thus cannot contribute to any transport. The quantum critical 
point at which the long range order appears has many features that distinguish it from a 
conventional quantum critical point in a translationally invariant systems, most notably it 
is characterized by a wide distribution of the order parameter in a realistic system and the 
appearance of a new intermediate phase in which only low energy local excitations have 
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infinitely long lifetime while high energy excitations can decay. 

In Section [TTl we analyze the experimental data on superconductor-insulator (SI) tran- 
sitions in disordered films of InO, TiN and Be and argue that this transition is driven by 
the competition of disorder and superconductivity with very little effect of the Coulomb 
repulsion. This will allow us to formulate a theoretical model for this quantum transition 
which is nothing but the model introduced in a seminal work by Ma and Lee^. In Section 
mil we develop the formalism to study the formation of the order parameter in this model 
at zero temperatures. In this formalism the appearance of a wide distribution of the order 
parameter shows up as replica symmetry breaking. Our theory provides the justification for 
the qualitative idea of the importance of rare sites characterized by a very large suscepti- 
bility, an idea that was first suggested by by Ma, Halperin and Le^ and observed in the 
solution of similar one-dimensional modelsS^S. In Section IIVI we study the properties of the 
insulating state. We first determine the level width (decay rate) at zero temperature and 
then extend the analysis to low temperatures. We find two phases in the resulting insulator: 
an intermediate phase where only excitations of large enough energy can decay and a third 
phase, in the strong disorder regime, characterized by the infinite lifetime of all excitations, 
similar to the one proposed in^. In Section |V] we discuss the effect of a magnetic field on the 
phase diagram within our model of the SI transition. The main conclusion of this Section is 
that the effects of the frustration induced by magnetic field are small; this can be understood 
as a consequence of the strong inhomogeneity of the order parameter in the vicinity of the 
transition. In Section |VT] we discuss the direct implications of the theory for the experiments 
and propose numerical simulations that should test the applicability of the theory to realistic 
models. Section |VTT] gives conclusions. 

A first quantitative study of the phase diagram of the Ma-Lee model has appeared recently 
in^. The present paper provides a much more detailed derivation, and studies in detail the 
behavior of the order parameter and level width. Similar qualitative conclusions on the 
reelvance of the Ma-Lee model and the phase diagram were reached from phenomenological 
considerations in the recent paper-. 
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II. SUPERCONDUCTOR-INSULATOR TRANSITION: DATA AND THE 
MODEL. 



A. Experimental results 

We begin with the analysis of the data on strongly disordered superconducting films. 
Our goal is to formulate the simplest model that captures the essential physics of these 
systems. Strongly disordered films of InO, TiN or Be display a zero field transition from 
superconductor to insulator when their resistivity in normal state exceeds a value of the order 
of the resistance quantum Rq = 6.5kfl.--^ The films in the superconducting state close 
to the transition become insulating when subject to a magnetic field.—"— The transition 
driven by magnetic field display a quantum critical point behavior: the resistance of the 
films at fields B < decreases with temperature decrease while the resistance of the films 
at B > Be increases with temperature decrease. Very important information is provided 
by tunneling spectroscopy of such superconducting films in the vicinity of the quantum 
critical pointr^"— . These data show a well defined gap at all points whilst the coherence 
peaks expected for a BCS superconductor appear at some locations and do not appear at 
others; a similar phenomenon was reported for high Tc oxides^"—. The absence of coherence 
peaks combined with intact superconducting gap in a single electron tunneling experiment 
implies that the disorder does not destroy local Cooper pairing of electrons but prevents 
formation of the coherent state of these pairs. This allows to exclude a fermionic mechanism 
of superconductivity suppression in these materials: in such an alternative scenario, the 
main role of the disorder would be to enhance the Coulomb interaction that competes with 
phonon attraction. This would lead to the suppression of the transition temperature and the 
superconducting gap and to the eventual disappearance of both the superconductivity and 
the gap. In contrast, the InO and TiN data demonstrate the presence of a large one-particle 
gap in the absence of global coherence and a smooth crossover between superconducting and 
insulating gaps as disorder is increased (see^ for more detailed discussion). 

In the absence of single-electron excitations, the supercondutor-insulator (SI) transition 
might happen either because the Coulomb interaction between the pairs prevents the for- 
mation of the condensate, or because the disorder of Cooper-pair energies prevents their 
coherent motion from one site to another. The first scenario might be realized in the granu- 
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lar materials in which superconductivity remains basically intact in each grain. The grains 
are coupled to each other by Josephson couplings which compete with the Coulomb energy 
that changes when a pair moves from one grain to another. Exactly the same physics is 
realized in Josephson junction arrays which were extensively studied some years ago.— "^^^^ In 
particular, it was established that in the presence of a magnetic field the Josephson arrays 
display a temperature-independent resistance that varies by many orders of magnitude as 
a function of a weak magnetic field. In other words, in these arrays the SI transition does 
not happen directly, instead there is a wide region of intermediate 'normal' phase. This 
behavior is very surprising given the absence of single electron excitations in these arrays. 
It is probably due to the formation of a Cooper-pairs glass, similar to the electron glass; in 
this regime collective modes provide the dissipation mechanism. Although the theoretical 
picture of this phase is not clear, the observed experimental behavior of these systems is 
in a striking contrast with the behavior of disordered films which show a direct transition 
between superconductor and insulator. 

This leaves the only possible mechanism for the superconductor-insulator transition in 
disordered films: the competition between pair hopping and random pair energies on different 
sites, as suggested 25 years ago in a seminal paper of Ma and Lee^ (see also^^i^— ). As we 
show in the present paper, the solution of this model reproduces correctly the most important 
features of the data: direct SI transition, activated behavior close to the quantum critical 
point in the insulating phase, strong dependence of the activation energy near the quantum 
critical point and huge order parameter variations from site to site in the superconducting 
phase. 

Recent experimental data indicate the possibility of a faster than activated temperature 
dependence of the resistance in the insulating phase (characterized by growing T In i? at low 
T). This behavior is very unusual for disordered electron systems which typically display 
either activated, Efros-Shklovskii or Mott behavior characterized by Ini? oc T~" with a < 1. 
It can be understood if the pair excitations remain localized in space at low energies but are 
delocalized at high energies with the temperature dependent mobility edge that separates 
them. 

For the microscopic justification of this mechanism one needs to find the reason why 
Coulomb repulsion does not play a role in the superconductor - insulator transition in some 
disordered films. Phenomenologically, it is known^° that the dielectric constant in these 
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materials remains very large, k > 30 deep in the insulating phase of InOx. A large value of 
the dielectric constant between low energy electrons close to the Fermi surface allows one to 
neglect the effect of Coulomb interaction on the pairing even for relatively strong disorder 
that results in wave functions localization. Because the host matrix dielectric constant can 
only increase with additional carriers, the Coulomb interaction between conducting electrons 
close to the Fermi surface is strongly suppressed in this materials. The microscopic reason 
for that might be a very strong energy dependence of the density of states (in particular, 
it was computed for a typical small cluster of TiN^^), this implies that although density of 
states exactly at the Fermi level is small and the states there get localized, there are many 
single electron states within lOmeV range that provide large screening of the bare Coulomb 
interaction. 

A detailed analysis of the competition between Cooper pairing and wave-function local- 
ization is provided by a recent paper—. In particular, it shows that global superconductivity 
can survive in a range of relatively strong disorder. Within this scenario the localized elec- 
tron pairs are formed at high temperatures when two electrons bind together (with typical 
binding energy Ap) on a weakly localized state. At lower temperatures, Tc < Ap, coherent 
hopping of these pairs leads to the long-range coherence and to the formation of a super- 
conducting condensate. Direct experimental confirmation of this scenario is provided by the 
scanning tunneling microscopy dataM"— . At even stronger disorder long-range coherence is 
never formed, the resulting state is insulating, characterized by a large single-particle gap. 

The paper— studied the properties of the superconducting and insulating phases far from 
the transition. In this work we study the properties of these phases in the vicinity of the 
transition in which the characteristic energy scales are much smaller than the single particle 
gap, Tc <^ Ap. In this study we shall focus on these low energy scales and neglect the contri- 
bution from single-electron occupied localized orbitals. This allows us to describe the physics 
entirely in terms of Anderson pseudospins^^. The interaction between the pseudospins is due 
to the matrix elements of the Cooper attraction between original electrons. As discussed 
in the work— these matrix elements are modified by the fractal nature of the electron wave 
functions which leads to correlations and to a smooth dependence of on disorder strength. 
In the present paper we shall ignore the effect of fractality that should not affect the qual- 
itative properties in the vicinity of SI transition. In the absence of such correlations the 
average interaction between pseudospins does not depend on disorder. Thus, when compar- 
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ing with experimental data we shall assume that changing of the disorder translates only 
into changing the number of "interacting neighbours" of a given pseudospin. 

B. The model: localized pairs and pseudospin representation. 

In the absence of Coulomb repulsion the only interaction between electrons is the attrac- 
tion that leads to their pairing on the localized single electron states.^ The strength of this 
pairing is inversely proportional to the volume of the state. Because of the fractal nature 
of the single electron wave functions the volume of a typical localized state is small, this 
makes the pairing energy large. Thus, the low energy degrees of freedom in this system are 
pairs that can hop from one site to another. This physics is described by a Hamiltonian 
of disordered hard-core bosons, or, equivalently, pseudospin operators (originally introduced 
by Anderson^^ for the pure system, and later by Ma and Lee^ for the disordered system): 

HxY = -2 5^6< - $^M,,(4s- + srs+) (1) 

where Sj = |crj are spin-| operators and the sum Yl(ij) goes over all different pairs of 
neighbors The state with = ±| corresponds to a local level occupied or unoccupied 
by a Cooper pair; ^^'s are occupation energies for each site, which are quenched random 
variables drawn from a probability distribution p{C,)- Hereafter we assume a box distribution 
P(0 = 2W^^(^ ~ I^D' ^^^^ non-interacting density of states is u = 1/2W. The 
important feature of this distribution is that it is constant near to ^ = 0; the value of u just 
sets the scale of energies, and we shall choose u = 1. The matrix elements Mjj describe the 
hopping amplitudes of Cooper pairs. These hopping amplitudes couple a typical local level 
to a large number of neighbors, Z ^ 1. We shall assume that each site is coupled to Z 
neighbors with Mij = 2g/{Z — 1). Another closely related problem corresponds to the Ising 
ferromagnet in a random transverse field; the model is defined by the Hamiltonian 

Hi = -Y.^^-- - ^Y.^--] (2) 

For brevity we shall refer to this model as Ising model below. We shall mainly study the XY 
problem ([T]) but most of our conclusions also hold for the case (E]). As will be clear below, 
in the leading approximation the two models ([1]) and ([2]) lead to identical cavity equations 
for both the order parameter and the relaxation rate. 
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The Hamiltonian ([T]) describes the superconductor-insulating transition as a ferromag- 
netic spin i system with random transverse fields, as proposed originally in^i^. In this 
language the superconducting phase corresponds to the existence of a spontaneous magne- 
tization in the x — y plane. 



III. FORMATION OF A LONG-RANGE SUPERCONDUCTING ORDER. 
A. A simple mean-field approach 

The most obvious approach to study the Hamiltonians (fT][2|) is to use the mean-field 
approach, in which Hxy is replaced by Hmf = ~ Ba^) where B is determined 

self-consistently from the equation B = {g/{Z — 1)) J^ji'^j)- The right-hand side of this 
equation contains a sum of Z ^ 1 terms, so it might be expected to become site-independent 
in the large Z limit. Although we shall challenge below the assumptions and conclusions of 
this approach, it is instructive to summarize its main results. The mean field Hamiltonian 
on site i has two eigenstates with energies — Tj a/ + -B^ where Ti = ±1. At temperature 
T = one obtains the self consistent equation for B, which gives either B = or 



.-,f^m '-''^^-\ (3) 

At zero temperature this equation has a solution for any g, so the ground state is char- 
acterized by a non-zero field, 5 > 0. This corresponds to a spontaneous order parameter 
(in spin language, this is a non-zero magnetization) in the x direction, so the system is a 
superconductor (or a ferromagnet in case of the Ising model). At finite temperature there 
is a transition from an insulating to a superconducting phase, when (j3]) starts to have a 
solution. Assuming a constant density of states, p(^) = \0{1 — |^|), we find that the critical 
temperature is determined by: 

1 ir'^^U^ ^ T.^l£^exp(-1A, (4) 

where C ~ 0.577.. is the Euler number. As we show below, while this mean field prediction 
is correct at Z = oo, it is qualitatively wrong for a system with a finite connectivity (even 
if Z is large) in the limit of very small g. 
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B. Cavity mapping 



We now turn to a more refined mean-field discussion, valid for finite Z S> 1, which is the 
basis for our main results. We shall employ the cavity method that has been developed in 
the classical case to study frustrated systems on a Bethe lattice, a fixed connectivity random 
graph of connectivity Z, seer^. Its generalization to quantum problems without disorder has 
been studied iiM-^, using the Suzuki- Trotter formalism. We shall use here an approximation 
of the quantum cavity method which does not use the Suzuki- Trotter formalism, and allows 
to study analytically the quantum disordered systems. As we shall see below (see discussion 
below (pTl) and section [Til Gp . this method takes into account the leading corrections to the 
naive mean field at large-but-finite Z, which turn out to be of order 1/lnZ . It neglects the 
corrections that are small in 1/Z. Appendix A explains the relation of our approach to the 
full Suzuki- Trotter quantum cavity method. 

In the simplified cavity method one emulates the effects of spin environment by the static 
field acting on it. One thus studies the properties of a spin j in the cavity graph where 
one of its neighbors has been deleted, assuming that the K = Z — 1 remaining neighbors 
are uncorrected. The system of spin j and its K neighbors is thus described by the local 
Hamiltonian 

K 
k=l 

where B^. is the local "cavity" field on spin k due to the rest of the spins (in absence of 
j). Here we have chosen the x direction as the transverse direction where the spontaneous 
ordering takes place. By solving the problem of Z spins in (j5]), one can compute the induced 
magnetization of j, which is by definition equal to ^^f^^^ tanh/3^^| + Bj. One thus 
gets a mapping allowing to compute the new cavity field Bj in terms of the K fields B^ on 
the neighboring spins. This mapping induces a self-consistent equation for the distribution 
of the B fields-^. Notice that the use of the mean-field approximation in order to study the 
cavity mapping is very different from its naive use in ([2]). 

One can make one more approximation which simplifies the cavity mapping and makes 
possible an analytical study. Similarly to the mean-field approach discussed above, one can 
replace dynamical spin variables by their averages. This approximation neglects some 
terms of order of l/K; below we refer to it as the cavity- mean- field approximation. We 
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shall first study this approximation here. Later on, in section IIII Gt we shall compare its 
results with the ones found by the full cavity mapping. As we show in this section the cavity- 
mean-field approximation becomes very accurate even for moderately small K. The physical 
reason for this is that the main effect missed by the cavity-mean-field approximation is the 
level repulsion induced by the interaction between sites j and k. However this repulsion 
happens only in the rare cases when is close to ^j, and it has a significant effect on the 
resulting value of the local order parameter only when both of them are close to 0. Because 
this happens very rarely, the cavity-mean-field approximation turns out to be very good for 
this problem. 

Formally, the cavity mean field approximation amounts to approximating the cavity 
Hamiltonian ([5]) acting on spin j by 

K 

Tjcav—MF c 'X 9 \ ^ / „x\ fa\ 

= -^j-^i -^j-;^Z^(^fc) (6) 

k=l 

This implies that Bj = J2k=i('^k)^ giving the recursion equation relating the B fields: 

^.■ = |E7J=ft-h/3^^|7^. (7) 

k=i V^fe + ^k 

This recursion induces a self-consistent equation for the distribution P{B): 
P{B) = j dB,... dBKP{B,) . . . P{Bk)5 {^-j^jz tanh + . (8) 

This equation always allows the trivial solution P{B) = 6{B) corresponding to an insulator. 
The superconducting transition is signaled by the appearance of another solution character- 
ized by a non-trivial distribution function P{B). It turns out that this transition and the 
properties of the phases in its vicinity can be best studied using methods developed in the 
statistical physics of random systems. 

C. The directed polymer problem and replica symmetry breaking 

In order to understand the properties of the mapping ([7]), let us imagine that we iterate 
it L ^ 1 times on a Bethe lattice. For L finite and N oo the corresponding graph is just 
a rooted tree with branching factor K at each node and depth L (see Fig. [1]). The field Bq 
at the root is a function of the fields on the boundary. In order to see whether there 
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FIG. 1: Iteration of the cavity equations three times (here with K = 2). The field Bq on the root is 
proportional to the infinitesimal field B added on the boundary. The corresponding susceptibility, 
given in ([9]), is equal to the partition function of the DP problem, obtained by summing over all 
the paths from the root to the boundary (such as the one shown by the dashed line). In this sum, 
each path has a weight equal to the product of terms associated with each edge that it visits. 

is spontaneous ordering, we study the value of Bq in linear response to infinitesimal fields 
Bi = B 1 on the boundary spins. This is given by 

g tanh(/3^fc) ] 
a J ■ 

where the sum is over all paths going from the root to the boundary, and the product 
]^,„gp is over all edges along the path P. The response S is nothing but the partition 
function for a directed polymer (DP) on a tree, where the energy of each edge is e~'^* = 
{g / K){tsinh.{f3^k) /^k) and the temperature has been set equal to one. The general method 
for the solution of such problems has been developed by Derrida and Spohn^. The solution 
can be expressed in terms of the convex function 

i, f (^£f^ip)] . (10) 

Let us denote by a; = m the value of x where this function is minimal. In the large L limit, 
there exist two regimes for the DP problem: 

• "RS-phase": If m > 1, then (l/L)lnS = /(I) + \'n.{g/K). In this case the supercon- 
ducting phase appears at a value gc = Ke~-^^^^ which is the same result as found by 
the naive mean field approach in (jl]). 
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• "RSB-phase": If m < 1, then InS = f {m) + \n{g / K) . The ordered phase appears 

at gc = Ke~f^'^'> which is larger than the naive estimate (jl]). 

The method for obtaining these results developed by Derrida and Spohrt^^ employed a map- 
ping to the traveling wave equations. An alternative way, which we now summarize, is to use 
the replica method, similar to the one used ivM. It gives a compact solution and helps to un- 
derstand the physical meaning of the DP phase transition. The names of the two phases used 
above ("RS" stands for replica symmetric" and "RSB" stands for replica-symmetry broken") 
originate from this analysis. 

The DP partition function S, defined in ([9]), depends on the random quenched variables 
^n- One expects that the free-energy in a short-range-interaction problem is a self- averaging 
quantity, so the value of log H for a typical sample is obtained from the quenched average of 



log S over these random variables, denoted by In S . In the replica method one computes it 
by writing 



InH = lim(S" - 1) 

n->0 



a=l..n \ Pa kdPa 

The average of S" is obtained by a sum over n paths, 



T=^i[U'-^^y, (12) 

Pl,...,P„ k ^ J 

where the weight of each edge k in the tree depends on the number of paths which go 
through this edge. 

The RS solution assumes that the leading contribution to f|T2|) comes from non- 
overlapping independent paths (r^ = 1). This gives 



^|tanh(/30 



exp(Ln[log((7/ir) + /(l)]) = (S)" . (13) 



The RSB solution assumes that the leading contribution to (fT2|) comes from patterns of n 
paths which consist of n/x groups of x identical paths, where the various groups go through 
distinct edges. This gives: 



- — - L ( g tanh(/3^fc; _ 
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ln(|)+/(x) (14) 



where f{x) is the function introduced in (fTOi) . In the rephca hmit n — )■ 0, the parameter 
X should belong to the interval [0,1]. Minimizing the function f{x) over x G [0,1] (the 
fact that one should minimize /, and not maximize it, is a well-known aspect of the replica 
method^^), one gets the phase diagram described above: there exists a critical value of 
the inverse temperature /3 = Prsb such that, for {3 < Prsbi the function f{x) is minimal 
at the boundary, x = 1, of the interval [0,1]; the DP is in its replica symmetric phase. 
For /3 > Prsb, the function f{x) has a minimum inside the interval (0, 1) at some value 
X = m < 1, this corresponds to the spontaneous breakdown of the replica symmetry in the 
DP problem. 

These two regimes of the DP problem are qualitatively very different. In the RS regime 
the measure on paths defined in iQ is more or less evenly distributed among all paths. On 
the contrary, the RSB regime is a glass phase where the measure condenses onto a small 
number of paths. An order parameter which distinguishes between these two phases is the 
participation ratio Y = "^pWp, where wp is the relative weight of path P in the measure 
([9]). It is easy to see that Y = in the RS phase. In the RSB phase, the value of Y is finite 
and non self-averaging (it depends on the realization of the ^'s), and its average is given by 
1 — m. This glass transition, and the nature of the RSB glass phase, are identical to the 
ones found in the random energy model^^i^. 

D. Phase diagram 

The results on the DP problem described in the previous section allow to derive the 
phase diagram of disordered superconductors. The state of the model is determined by the 
three parameters: the coupling g, the cavity connectivity K = Z — 1, and the temperature 
T = 1//3. The phase transition between insulator and superconductor is a surface in this 
three-dimensional parameter space. Depending on the purpose, it can be useful to view 
it in different directions. We shall define the phase transition as gc{K,T), or Kc{g,T), or 
Tc{K, g) = l/Pc{K, g). In the zero-temperature limit we shall speak of gdK) = gc{K, T = 0) 
andK,{g) = K,{g,T = 0). 

The DP partition function S gives the susceptibility, measured on the root site, to a small 
field on the boundary. When this susceptibility diverges, there is spontaneous generation of 
non-zero B fields, i.e. the systems is in the superconducting phase. The superconducting 
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phase transition point is given by: 

• If /3 < I3rsb (RS phase of the DP), then gc = Ke^f^^\ This coincides with the result 
of the naive mean field analysis. 

• If /3 > Prsb (RSB phase of the DP), the phase transition is obtained by solving the 
two equations f{m) + In = and f'{rn) = 0. 

At zero temperature the function f{x) can be computed analytically: 

fix) = i In (-^] (15) 



X \1 — X ^ 

leading to two equations for m and Qc which have the explicit solution 

m = 1 — egc (17) 

The value of Qc determined by the solution of f lT6|) gives the location of the zero temperature 
quantum phase transition gc{K) between the insulating and the superconducting phase. 
Note that in the framework of the recursion equations ([7]) this is an exact result. It predicts 
a finite value of Qc for all i^' > 1. In particular, when = 1 it reproduces correctly the 
result Qc = 1/e known for the one- dimensional chain^. This is in contrast with the naive 
mean field equation (jl]) which wrongly predicts Qc = for all K. For the comparison with 
the experiment it will be useful to describe the phase transition point as a function of the 
effective number of cavity neighbors K for a given a value g of the hopping (see section Hi Al) ; 



in terms of this parameter the system becomes superconducting at K > Kc{g,T). The zero 
temperature transition point is given by explicit equation: 

K,{g) = ge'/^'^\ g<l/e. (18) 

We now discuss the finite temperature phase transition. At sufficiently high temperatures, 
T > Trsb = ^/f^RSB, the transition line Tc{K,g) is independent of K and is given by the 
naive mean field result (jl]): Tc = (4e*"/7r)e^^/^. The RSB solution appears at a value gnss 
such that Tc{K, gnsB) = Trsb- At this g the derivative /'(I) = 0, which gives two equations 
that determine the position of this point on {g, T) or {K, T) plane: 

(7 tanh/3e = 1 (19) 

g [ ^In ( -l^tanh/Jn tanh/3^ = (20) 



e 
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The result is again best expressed by inverting gjissiK) in order to get the critical value 
KRSB{g) below which the naive mean field solution breaks down. The asymptotic of the 
solution of the equations at 5^ <^ 1, K ^ lis, 

ir^5B(^)^^gl/(2,) (21) 

The equations for the critical values of the number of neighbors (118f2ip show that, at large K, 
the solution deviates from the naive mean field solution of section UlI Al at g < l/(elni^). 
This demonstrates that our approximation takes into account 1/lni^ corrections, as an- 
nounced above. Note that at weak coupling K^^^{g) is exponentially larger than Kc{g), so 
that there is a broad range of K where RSB effects are important. 

An important unexpected qualitative conclusion from (12T1) is that the effective number 
of neighbors of a given site, Zeg, defined as the number of neighbors having local ener- 
gies within the energy stripe ±Tc, is typically much smaller than unity. Indeed, at the 
temperature Trsb, this number is given by 

Zeff = K^'^T^o = -e^'ge-'l^' « 1, (22) 

TT 

it decreases with K roughly as Z^s ocl/K when K ^ 1. Although Z^ff is small, it leads to 
a non-zero T^. and to RSB effects. This is another signature of the fact that the transition 
is governed by rare sites. 

The phase diagram for = 4 is shown in Fig. [2l At any temperature, there is a finite 
critical value of the coupling, gc{K, T), separating an "ordered", superconducting, phase with 
spontaneous order parameter aX g > gc{K,T) from a "disordered" normal phase with zero 
order parameter at g < gc{K,T). Contrary to the naive MF prediction, gc{K,T) remains 
non-zero at T = 0. 

One of the important characteristic properties of the quantum critical point is the behavior 
of Tc{K,g) near to the critical point Kc{g) = ge^/^^^^ where it drops to zero. In order to 
study this behavior, we look for the solution of the two equations /(x) + \og{g/ K) = and 
f'{x) = 0. The explicit form of these equations is 

<if^^lnf^tauh,3?Uo (24) 
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FIG. 2: Main panel: phase diagram in plane {g, T) for K = A. The full lines show the critical 
temperature as function of g. The low temperature phase is superconducting, the high temperature 
phase is insulating. The top curve is the naive mean-field prediction which gives the correct result 
above Trsb = 0.0207. The bottom curve is the result of the analysis on the Bethe lattice described 
in the text that includes the RSB effects in the DP problem, which occur at temperatures T < Tusb- 
The insert shows the critical temperature as function of K for g = 0.129, a value which roughly 
corresponds to the experimental situation in disordered InO films (see section rvT|) . For this value 
of g, the replica symmetric solution gives a i^-independent transition temperature = 0.001. 
This prediction of the replica symmetric theory is correct for K > K^^^ ~ 6. For smaller K the 
transition temperature starts to drop, the quantum critical point corresponds to Kc — 2.2. Notice 
that in a rather wide regime the replica symmetry is broken but the effect of the breaking on the 
transition temperature is small. 



We assume that K = Kc{g){l + 6) with 5^1, and we introduce the quantity y = 1 — x, 
which is expected to be close to its critical value i/c = eg. The second equation ( l24l) leads 
to the relation y{6) = y^l — yc^)- Next, in order to determine Tc{6) we use the stationarity 
of the first equation with respect to the simultaneous variations of T, K and y, and find 



with 



T,{K) = i}iy,) ( ^ - 1 



C{y) 



L(i-y)c(2/), 



(25) 



dt 



t 



tanht 



cosh t 



The result ( l25l) is valid as long as Tc{K) is much less than the mean-field transition temper- 
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ature T^o ~ e and K — ^ K^. In terms of K — the former condition is 



For physically relevant systems g ^ 0.02, so the condition f l26|) determines the regime of i^' 
for which the result ( !25ll holds. This analysis shows that the transition temperature Tc{K) 
goes down very slowly when K decrease below the K^^^{g) value, until it reaches the close 
vicinity of Kc ^ K^^b ^j-^g^-g it drops sharply to zero according to Eq. fl25|) . The numerical 
solution of Eqs. (12311241) for Tc{K,g) at g = 0.129 is presented in Fig|2] (insert). 

The replica symmetry breaking at i^' < K^^^, which is at the origin of the failure of 
the naive mean field analysis, has important physical consequences. Physically it signifies 
the absence of self averaging which is due to the importance of very rare sites with small 
C,j that easily polarize in the x-direction, becoming thus nucleation centers for the ordered 
phase. The qualitative importance of these rare sites and resulting absence of self averaging 
was first discussed in~. Our analysis gives the mathematical formalism to evaluate these 
effects and their consequences within the well-controlled Bethe approximation, such as the 
exact position of the phase transition point. The role of rare sites with small C,j can be 
quantitatively characterized by the condensation of the measure in the DP problem. Let us 
imagine adding a small field in the x direction on one site j, and let us look at its effect on 
the spins k at distance L from j; this gives the picture of the correlations in the many-body 
wave-function. In the RSB phase, the order develops only along a small number of paths; 
this fact is completely missed by the naive mean field approach. 

Even at > K^^^ the self averaging is not fully restored. In particular, in a wide 
range of K the simple mean-field approximation is not applicable for the calculation of 
higher moments of the P{B) distribution. The upper border of this region is given by the 
value K* corresponding to the divergence of infinitely high moments (divergence of at 
n — oo). The condition for this divergence is given by the inequality g/KT^ > 1. With Tc 
determined by the mean-field approximation, we get the critical value K* that corresponds 
to the divergence of large moments: 

^ tanh z = l ^K* = -e-'^ge^l^ (27) 
where C is the Euler number. The full self- averaging is expected only at K > K* . 
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E. Distribution function of the local order parameters and the phase diagram. 



In this section we derive the equation for the distribution function P{B) of local fields 
close to the transition point, T ^ Tc{g,K), using the cavity mapping ([7]) as the starting 
point. We will use the Laplace transform of this distribution, V{s) = dBP{B)e^^^ . 

We start from the linearized version of Eq.([7]), 

k 

which is adequate to describe the phase transition region where the fields are small. The 
Laplace transform satisfies the equation: 



^^ = |Er*'^''^^^^') ' ^^^^ 



Vis) 



\ ( ^tanh^r^"" 



(29) 



Consider the first terms of its expansion at small s, and assume that this expansion behaves 
as V{s) = 1 — As^ . If X = 1, the mean j dBP{B)B is finite. We will see that this is the 
situation in the RS phase. The situation x < 1 occurs in the RSB phase, and corresponds 
to a distribution P{B) which decays at large B like B~^~^, with a diverging mean. 

Let us first assume that x = 1. Then A is just the mean of B, i.e. it coincides with the 
usual superconducting order parameter. Plugging V{s) = 1 — As into (!29ll . the resolvability 
condition for the linear equation in A leads immediately to the standard mean-field equation 
for Tc, see (jlj). We learned previously from the replica analysis that this equation is valid 
only in the RS regime: it is not correct when K < Kc{g). Therefore, RSB corresponds 
to the divergence of the mean order parameter. To study the RSB phase, we assume that 
V{s) = 1 — As^ with X < 1. From (!29|) . we obtain the condition 

(30) 

which can also be written by using the function / introduced in the replica analysis as 
/(x) + \og{g/K) = 0. Therefore, for any x < 1, there is a non-trivial solution g = Ke~f^^\ 
In this situation, it is reasonable to assume that the critical value of g is the largest one 
among all these values, which means that it is obtained by finding the minimum of /(x). 
We shall prove this assumption in section IIIIII where we study the spatial evolution of 
the distribution function. This is precisely the result obtained with the replica solution in 
section [Tll Dl of the DP problem, with x = m defined in Eq. (|T7|l . a solution which coincides 
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with the result obtained by travelhng waves analysis^^ when apphed to this problem. The 
corresponding distribution function P{B) is characterized by a power-law behavior at large 
B ^ Bq, where Bq is of the order of a typical value of B: 

PiB) = ^ (31) 
This power-law tail at large B translates into the behavior of the Laplace transform: 

V{s) ^ 1 - (sfio)™ at sBo < 1 (32) 



F. Scaling of the order parameter close to the transition 

The expansion of the mean field equation ([3]) would give an order parameter which scales 
as i? ~ V9 ~ 9c close to the transition. This is the usual mean-field scaling. It is modified by 
the strong fluctuations of local ordering flelds present in our problem. As we discuss in section 
IIIIF ll this modiflcation is present even in the some range of parameters in replica symmetric 
phase not too far from RSB transition and, of course, in the whole low temperature regime 
where RSB happens. We shall analyze these two cases successively. 



1. Replica symmetric regime at large K. 

The nonlinear recursion relation ([7]) leads to the equation for the average value 

{B)=gl dBP{B) £ di ^^_g__ tanh(/3v/g^) (33) 

At T > TjisB, or K > K^^^{g) defined in Eq. (l2T|) . the transition temperature is determined 
by Eq.(jl]). We assume now that, just below the transition line Tc{K,g), the distribution 
P{B) has a scaling form P{B) = Bq^^{B/Bo), with Bq vanishing at T -> T~. 

We begin by repeating the usual mean field arguments. Assuming that $(x) decays faster 
than so that the integral dx^{x)x^ < oo, we can expand 

tanh{(3^/WT^) tanh(/30 5^ / f3^^ 



etanh(/30- 



We use the equation (jlj) for the critical point to reduce fl33|l to the form 



Bq J rfxx$(x) (^^-l^ =Bl j dx^{x)x^ di^ {^i&n\i{(5Cj 



cosh'(/30 



(34) 
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which leads to the usual mean-field scaling Bq ~ y/g--^- 

Clearly, the crucial assumption used in this solution is that the third moment of P{B) (i.e. 
J ^{x)x^dx) is finite. It breaks down if P{B) decays at large B like B^^^"" with 1 < a < 3. 
Notice that in this regime the mean value of B is finite, therefore it belongs to the replica 
symmetric phase from the point of view of the equations giving the critical temperature 
discussed in section UlI D[ The behavior of P{B) at large B can be studied by the Laplace 
transform method developed in section 1111 E[ More precisely, we assume that the Laplace 
transform at small s has a form l — Bs — As"" with 1 < a < 3, insert it in the general recursion 
equation fl29l) and solve the resulting equations for the coefficient A and the exponent a. The 
computation shows that the exponent a is the solution a > 1 of the equation /(a) = /(I). 
The mean field scaling holds if and only if this solution satisfies a > 3. This happens when 

S1 1/2 

(35) 



K>K. = 



c 



The result 1 — Bs — As"^ for Laplace transform implies that in the regime K^^^ < K < 
the distribution P{B) on the critical line T = Tc{g) decays as CBq/B^^"- with a < 3 and 
C ~ 1. In this case the equation ( l34|) for the order parameter is replaced by 



^ ' 'Jo Jo \ ^Blx^ + e ^ J 



(36) 



This implies the anomalous scaling of the order parameter near the critical line: 

Bo{g) ^ (g - gcY^^--'^ (37) 

Note that when K — the exponent a — )■ 3 so that Bo{g) dependence reduces to the usual 
square-root singularity. At the opposite end of the interval K^^^ < K < K^, the exponent 
in Eq.( l37|) diverges because a ^ 1 as K }<[RSB _ 



2. Regime of broken replica symmetry. 

The RSB transition affects dramatically the scaling of the field in the ordered phase in 
the vicinity of the transition. As we have seen in subsection IIIIF 11 in order to derive this 
scaling we need to know the distribution of fields P{B) induced by the mapping ([7]). As 
discussed in section IIIIEI the expansion of its Laplace transform shows that exactly on the 
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transition line g = gc{K,T), in the RSB regime, the distribution function P{B) decays at 
large B as P{B) oc l/B^^'"^. Here m < 1 is the RSB parameter identified in our analysis of 
the directed polymer problem as the value of x where f{x) is minimal. 

For g larger than g^ but close to the transition, a small typical field, Bq, appears. At 
B ^ Bq the usual scaling arguments show that the distribution P{B) retains the same 
scaling form (!3T!) that it acquired at gc- As we show below this scaling law is cut off from 
above by Bmax — g/ K ^ Bq , beyond which P{B) decays rapidly. The power law behavior 
of P{B) in a wide range of fields allows one to find analytically the typical field Bq from the 
following reasoning. Because the integral J P{B)B^dB ~ \n{Bmax/ Bq) is logarithmic, the 
main contribution to the expectation value of -B™ comes equally from a very broad range of 
fields Bmax ^ B ^ Bq. Because of this wide range of fields contributing to the expectation 
value, the m'th power of ([7]) is dominated by the largest term in the sum, so that 

^r^E(f^-K,y]^j . (38) 

Similar arguments of dominance of the largest term justify the cutoff aX B c:^ g/K ^ Bq. 
Because each term in the sum ([7]) is cut off by Bmax = g/K and the sum is dominated by the 
largest term, the power-law behavior ( 13T|) is effectively cut off by B^ax, while at i? > Bmax 
it decreases exponentially. 

Averaging the approximate mapping for 5™ (|38l) we get the equation: 

To proceed further we write the equation for the transition line, gc = Ke~^^"^\ in the form 
similar to Eq. fl39|) : 

and subtract (gO]) from ([39]): 



]{B"') = K (^j^y P{B)B"'dB 



'tanh(/3v/52TF) V Aanh(/30\ 
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As will be clear below, the qualitative properties of the solution are the same at zero and finite 
temperatures. So, to simplify the formulas we focus on the vicinity of the zero-temperature 
quantum critical point Tc = 0, g = gdK). In this regime we can replace all tanh(/3M) 
functions by 1 and get 

9 



where 



7(x) 



dt 



fx (^2 + l)x/2 



In the large K regime which is of main interest to us, x — )■ 1, and 7(x) — i- (1 



X 



,-1 



(42) 



(43) 
. Using 



our Ansatz for the distribution function, valid in the broad range of B^ax ^ B ^ Bq, we 
find: 



1 - 



,„ *-(|)'"SWr;; (44) 

Here, and in the following, C denotes a numerical constant of order 1. Using the estimate 
for Bmax discussed above we get 

C7(m) 



Bo{9) 



~ e ^/'-''S'^) exp 



For e 



{g/g^)m-l_ 

— 1 ^ 1 we expand the exponent in Eq. fj45|) and find 

C7(m) 



Boig) 



~ e~i/(^^^) exp 



(45) 



(46) 



m{g - gc 

This gives the typical value Bq of the order parameter close to the transition; in contrast to 
quantum critical points in clean systems the behavior -Bo(fi') displays an essential singularity. 
This result becomes even more surprising when one compares Bq to the critical temperature 
Tcig) = C{g — gcY^^^^^ (which follows from Eq.f l25l) ) because it implies that the typical order 
parameter at T = is much smaller than the transition temperature at the same value of 

9- 

Eq. fH5|) gives the dependence of the order parameter as a function of the interaction 
constant at very low temperatures for g close to gc- Similar arguments give the dependence 
of the order parameter at (7 = gRSB as function of temperature in the vicinity of transition: 



BoiT) 



c 



(47) 



2gRSBHTc/T)\ 

Because this analytical derivation relies on our Ansatz for the distribution function, we 
have checked the scaling behavior (jl5]) numerically. We solve the self consistent equation IM) 
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for P{B) using the population dynamics algorithm developed in work—. In a nutshell, this 
method amounts to approximating P{B) by a population of fields, which is sampled by a 
Monte Carlo method iterating Eq.Q. In order to tame the effects of rare fluctuations, it is 
numerically convenient to put the system in a small external field, using the mapping f HH]) 
defined below in the Section IIII HI The results are reliable (as can be checked from the fact 
that they do not depend on the external field at small enough field) when g is not too close 
to Qc- The results for K = A (for which the zero-temperature critical coupling is Qc — 0.1 
and m ~ 0.73) are shown in Fig. [3] where we plotted log-B as a function of {g/g^^ — 1 
ioT K = 4 and its fit to a [{g/gc)"' - 1]"^ dependence: log 5 = A + EKg/gc)"" - 1]"^ with 
A = —4.0, B = 1.9. This value of B is close to the one expected from (l45ll because for this 
value of K one has 7(m) ^ 2.8. 

-10 
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-16 
-18 
-20 
-22 

0.10 0.15 0.20 0.25 0.30 0.35 

FIG. 3: The typical order parameter vanishes with an essential singularity at the superconductor- 
insulator transition. Upper curves show the logarithm of the typical fields (defined by In Btyp = 
InB]) obtained by population dynamics with small external fields: h = 10^^, 10^^, 10~^^ (from top 
to bottom). The lower curve shows the fit to l/{{g/gc)™' — 1) dependence expected from ()15]) . 

The computations above are based on the assumption that the distribution function 
retains its power-law dependence in the ordered phase in a wide range of fields. We have 
checked this assumption directly using population dynamics. The result is shown in Fig. HI 
One observes that P{B)B^^"^ develops a plateau that becomes wider on a logarithmic scale 
when one approaches the quantum critical point. 
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FIG. 4: Distribution function, B^~^'^P{B) for K = A and different values of {g — Qc)/ Qc shown next 
to each curve. As expected, B^'^"^P{B) is field-independent in a range of fields which is wide in 
logarithmic scale, and becomes wider as the transition is approached. 

G. Leading correction in 1/Z to the RSB solution. 

As we explained in section IIIIBt the full cavity Hamiltonian (|5]) can be approximately 
replaced by the simpler Hamiltonian (|6]) in which the dynamics variables were replaced by 
their averages (crj). This approximation ignores dynamic quantum correlations between the 
spin at site i and its neighbors at site j . These correlations might become important when 
the energies of these two sites are close, so that — g/ K . Such resonance conditions 
occur with a probability p ^ g/K which is small at large K. To check that these and other 
corrections are not relevant at modest values of K we have performed a population dynamics 
simulation of the full cavity mapping (|5]). At each step of this simulation we start with the 
population of spins characterized by energies and fields Bk. We randomly choose A^ 
sets of K spins from this population, add to each of them an additional spin with random 
energy ^ and diagonalize the corresponding Hamiltonian ([5]) to determine the effective value 
of the field acting on the additional spin. This gives a new set of spins with energies and 
fields -Bfc. A typical simulation involved 10'* spins and 10* steps which is sufficient to get the 
convergence of the typical field at {g — gc)/gc > 0.2. 

In Fig. [5]we compare the results of the full diagonalization with the simplified mapping for 
A' = 4. As one can see from these results, even at these modest value of K the results of the 
full cavity mapping are indistinguishable from the simplified mapping. Qualitatively it means 
that the important physics driving this transition is the appearance of wide distribution of 
fields, and this effect is not sensitive to quantum correlations neglected by the simphfied 
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FIG. 5: Comparison between the typical fields (measured by the harmonic mean) computed with 
the population dynamics of the full cavity Hamiltonian ([5]) and the simplified one ([6]). The line 
corresponds to the results of the simplified mapping, the dots to the mapping given by the full 
cavity Hamiltonian. 

mapping. 



H. Susceptibility in the disordered phase. 

We now study the approach to the quantum critical point coming from the insulating 
phase. In usual phase transitions the order parameter induced by a small external field 
diverges at the transition, corresponding to a divergent linear susceptibility. We shall see 
that the situation is very different here: there is no diverging susceptibility. 

We focus on the insulating phase at T = and g < Qc- We apply a very small uniform 
external field h to all sites. Repeating the same arguments as before we arrive at the cavity- 
mean-field mapping 

For very small h, in linear response, the induced fields are also small. We can thus neglect 
the non-linearity in this equation aX g < Qc- In the resulting linear equation we can rescale 
the variables Bi = bih and get 
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All properties of the solution are contained in the distribution function P{b) generated by 
the mapping fH9|) . In order to study this function, it is convenient to introduce the Laplace 
transform V{s) = f P{h) exp{—sb)db . It satisfies a self-consistent equation that can be 
derived directly from the mapping fj49|) : 



V{s) 



1 / 







(50) 



We assume the existence of V{s) and derive some of its properties. The behavior at small s 
is easily found: if we look for a solution in the form V{s) = 1 — as^, we get the equation for 
the exponent yu: 

ir^-V = 1 - 

This equation is equivalent to the equation for the RSB parameter m ai g = gc, see Eq. f|T6|) . 
Therefore fi = m. 

When s — 7- oo we look for a solution in the form V{s) = cs^exp(— as). Inserting it into 
the equation (!50|) we get, for s ^ K/g: 

f 1 .^y/--.X- ^^^^ 



1-9 K-1 \il-9) \9 

In order to prove the existence and stability of this solution we solved the equation f l50|) 
numerically. In order to check the analytical asymptotic we need to have the solution in 
a broad range of parameter s. For this reason we have transformed the equation fl50|) by 
introducing the variables s = e'', ^ = e~'^: 



cirexp(-r)p(log(-|)+C + r) (52) 

and solved the resulting equation by iteration for = 4 and g = gc = 0.0996. The result is 
displayed in Fig. [61 The low s asymptotic is indeed V{s) = 1 — as™" , with the correct value 
of m = 1 — 0.1 -e. The asymptotic behavior at large s is more tricky because the exponential 
dependence V{s) = cs*exp(— as) with parameters given in Eq. flST]) is expected to occur 
when gs/K ^ 1; this condition implies that in this regime V{s) ^ 10^^° which is difficult 
to access numerically. In the transient regime of K/g ^ s ^ 1 the numerical solution fits 
exponential dependence with somewhat different parameters, a = 1.2 and b = —1.7. 

The existence of a stable solution V{s) of the equation (!50|) implies a number of unusual 
properties of this T = quantum phase transition. A small external field h induces some 
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FIG. 6: Laplace transform of the distribution function produced by a small external field at g = gc 
for K = 4. Main panel shows the regime of small s, the inset shows the behavior at large s. The 
dashed lines show the asymptotic behavior at small and large s discussed in the text. 

non-zero fields Bi. Let us compute the average Bf. It can be expressed through the Laplace 
transform by 



The behavior V{s) = 1 — as"* implies that for x < m the average B^ ~ h^. In particular, the 
typical induced order parameter is of the order of the external applied field: exp(logi?) ~ h. 
But the moments B-^ with x > m, and in particular the mean B, are divergent at the level 
of the linear response to h. The non-linearity of the mapping neglected in (l49ll cuts off this 
divergence at B g/K. One can therefore expect a non-linear response B^ ~ {g / KY'^K^ 
when X > m. These results show that the response to an external field, computed at g < gc, 
has no singularity at the transition. This behavior is totally different from the one in usual 
phase transitions. We illustrate this by Fig. [7] which shows the average and typical fields 
induced by a small external field at g < gc- 

I. Spatial scale of inhomogeneities of the order parameter. 

Close to the transition the spatial scales beyond which the system is uniform become 
very large. In particular, at temperatures below that of replica symmetry breaking, the 
susceptibility is dominated by a single path, as discussed above, in Section IIII C[ implying 
that the system is essentially non-uniform at all length scales. The goal of this section is to 
compute the characteristic scales at which the system becomes uniform in the ordered state. 

The non- uniformity at short scales is related to the fact that close to the transition 
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FIG. 7: Average and typical fields induced by an external field at g < Qc and K = A. The lower set 
of curves shows the typical response Btyp/h = expilog B)/h for external fields h = 10"^, 10"'^, 10"^ 
and 10^^ (bottom to up). The upper set of curves shows the average response B^v/h = B/h for the 
same values of the external field. The curves show no singularity at the critical value of g indicated 
by arrow. As discussed in the text, the average response to the external field is controlled by the 
far tail of the distribution function and is non-linear, so the ratio Bav/h grows at /i — )• 0. 

the order parameter in the infinite system is power-law distributed in an exponentially wide 
range, from Bq to g/K. In contrast, at a given site the order parameter has some value which 
changes by a factor 0(1) in the vicinity of this site. Thus, at short scales the order parameter 
acquires values of the same order of magnitude, whilst the full distribution function is formed 
only at large scales. 

In order to describe this physics quantitatively, we write down the equations for the 
spatial evolution of the Laplace transform of the distribution function upon iteration on the 
Bethe lattice: 



J dbPniB)exp 



K 

9 sB ■ ■ 



(53) 



The crucial feature of the stationary solution of this equation is a power-law dependence of 
the distribution function in the exponentially wide range of s: 1 <^ s <^ ^/Bq. In this range 
we can neglect the non-linear (S^) term in the square root of (!53|) . the same is true for a 
general (non-stationary) solution of equation fl53|) in this parameter range. This allows to 
reduce the equation to the evolution of Laplace transforms: 



1 / 



(54) 



The stationary solution of this equation was discussed above, in Section IIII E[ Here 
we need to find the spatial scales (i.e. the number of iterations) at which this stationary 
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solution emerges for Pq{B) that corresponds to a particular value of the field, B = Bi, i.e. 
Viis) = exp{-sBi). 

The initial stages of evolution lead to the distribution function that has many features of 
the stationary solution given by Eqs. fl31|32p . V{s) = 1 — a(si?o)™', in particular it becomes 
close to unity in a broad range of s. The final spreading over the whole range 1 <^ s ^ l/i?o 
and thus the spatial scale at which the stationary solution is realized can be described by 
the linearized equation for = 1 — V{s): 



The evolution described by this equation approaches slowly the stable stationary solution 
0oo(s) = a{sBo)"^ found before. As we shall see below, this evolution is similar to a diffusion 
equation, so the total number of steps ("time") needed for this evolution is controlled by 
the final spreading of the distribution function (the smaller Bq, the more iterations it takes 
to reach the stationay solution). We can study this convergence towards (poo by assuming 
smooth deviations: 0ri(s) = ?/(ln s)0oo(s) , where y{() is a slow function of its variable. We 
get 



The integral over ^ in Eq. fl55|) is dominated by ^ ~ 1, whereas |ln((7/i^')| ^ l/{egc)- 
The assumption that y{() is a slow function on the scale of {egc)~^ allows us to expand 
j/n ( In + O in powers of In . Carrying this expansion up to the second order we get 



tice that, if this stationarity condition were not satisfied, one would get z/ 7^ 0, implying 
a non-zero drift term in the equation (!56l) . In this situation s stationary solution for the 
probability distribution would be impossible. This argument proves that the existence of 
a stationary solution of the recursion equation for the probability distribution of the order 
parameter implies that the value of m is fixed by the stationarity condition, as we found in 
the replica analysis. 

Altogether the integral equation (!55l) reduces to the equation of diffusive evolution: 





(55) 





Vn+liO -VniO = D 



■yniC) 



(57) 
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with a diffusion coefficient 



D^l^j-. (58) 



The longest relaxation "time" of this diffusive motion on the interval (0, In -^) is given by 



4 In^ 5, 



This relaxation "time" is actually equal to the correlation length of our problem. Close to 
the transition, Bq goes to zero as in Eg. (1461) . and therefore the length scale at which the 
system becomes essentially uniform diverges as N oc {g — gc 
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IV. WIDTH OF THE LEVELS IN THE INSULATING STATE. 

In the disordered phase the average value of the transverse field is zero. However the 
fluctuations of this field may be important. The main physical effect of these fluctuations is 
the broadening of the local levels that corresponded to = ±1 in the (? — )■ limit. We shall 
study this broadening both at T = and at T 7^ 0. Note that level broadening at T = 
is a rather complex phenomenon. It implies that a local excitation of spin i at frequency 
u = 2C,i decays. By energy conservation such a decay implies the excitation of some other 
spins. This cannot happen in a finite system because the energy of the spins are discrete 
and random. So the level broadening effect can only appear in an infinite system, where the 
excitation can propagate to an infinite number of other spins. We will show here that the 
broadening of levels in the insulating phase appears as a phase transition. 



A. Propagation of time-dependent perturbations, mobility edge. 

In order to study infinite systems consistently, we adopt an approach similar to the one 
developed above for the study of the transition into the ordered state. Namely, we consider 
an infinite Bethe lattice which is very weakly coupled to the environment at its boundary, 
study the effective level width at a distance L from the boundary and take the limit L — )■ 00. 
Thus we add to the Hamiltonian (j2]) the boundary term Henv = Ylii'^j'^ji't) + h.c.) where 
Xj{t) are dynamical fields, generated by the environment, characterized by a spectral function 
S{uj). In the leading order of the perturbation theory in g/K, the effect of these fields on 
the spin at a distance L from the boundary follows from the Fermi Golden rule. Imagine 
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that the environment of a spin j on the boundary induces a perturbation of frequency 00 in 



states of the spin corresponding to (Jq = ±1. This matrix element appears only in the L 
order of the perturbation theory in g and is equal to Xj{uj) / ^) / ~ 2|^fc|)], where 

the index k runs along the path connecting spin and the spin j at the boundary. Thus, 
in this approximation, the application of the golden rule to the relaxation rate of the spin 
gives, at zero temperature: 



where the sum runs over all paths connecting spin to the spins at the boundary at distance 
L. This equation is valid provided that all fractions inside the product remain small, the 
usual condition for the validity of perturbation theory. It should be modified when some of 
the fractions get large, but as we will see below these cases are so rare that these modifications 
are irrelevant. 

The study of the spontaneous emergence of a finite width can be done using the same 
directed polymer technique that we used in the previous section. Before we get into the 
details of this derivation it is useful to summarize the results that we shall obtain. When 
g < gd{oj, K) the levels have zero width, they are discrete. The curve g = gd{oJ, K) defines 
the spontaneous appearance of a finite width. (7^(0, K) is equal to the (zero-temperature) 
critical coupling where the system becomes superconductor: (7^(0, K) = gdK). For finite u, 
gdi^jJ, K) < gd{0, K); the minimal value of gd{oj, K) as function of uj occurs in the middle of 
the band, at w = 1. hi g < g*{K) = gd{u) = 1, K) the relaxation rate is zero for all states. 
This regime corresponds to the superinsulator introduced in^. For a given intermediate value 
of the coupling g*{K) < g < gc{K), the states in the middle of the band have a finite width. 
They are separated from the states of zero width by a critical energy Ud{g, K) (obtained by 
inverting the function gd{uj, K)) similar to the mobility edge of the non-interacting problem. 

We now derive these results using the mapping to directed-polymers. The partition 
function of the directed polymer on the tree is now: 



the form a^Xj{uj)e + h.c. This perturbation induces a matrix element between the two 




(60) 



En 



r 2g/K 



(61) 



p keP 
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The computation of log S is deduced from the evaluation of the function 

(62) 



1 / o \ 

K I 



f^{x) = - In . 

In the directed polymer terminology this problem turns out to be always in the low tem- 
perature phase corresponding to replica symmetry breaking, where the main contribution 
comes from a very small number of paths. This is due to the fact that fuj{x) always has 
a minimum at x = b < 1. As a consequence, one finds, using the same approach as in 



SectlmCl {l/L)\ogE = Ub) + 2\n{g/K). 

When oj = one has fuj=o{x) = 2f(2x), where / is given in ( fTOl) . Therefore, in the whole 



insulating regime g < Qc one gets (l/L)lnH < 0. Thus, the relaxation rate at very low 
frequencies decreases away from the boundary, in the bulk of the sample the width of levels 
is zero. 

Let us now study the case of non-zero frequencies. The critical line at which a finite 
level width appears in the bulk of the sample is given by the two equations ^ \x=b = and 
fu,{b) +2ln{g/K) = which read: 

9 



— ] / 1 Vt^t In — = In — (64) 

yKJ J, 1^-00/21^" \^-uj/2\ - ^ ^ 



These equations can be written explicitly: 

K'g^-' [(1 - uj/2y + {uj/2y] = z (65) 
iT"^^-" [(1 -tj/2)^ln(l -a;/2) + (a;/2)^ln(a;/2)] = 1-zln^ (66) 

where we have introduced the notation z = 1 — 2b. 

We begin by finding the region of the parameters {g, K) in which the system of equations 
f l63|M|) has a solution. Clearly, the solution of these equations expressed as a function of 
Kg{u) is symmetric under ^ — )■ 1 — ^ at fixed g. The energy u = 1 corresponds to the 



minimal value of Kg{l) = K*{g). For K < K*{g) the equations ( l63|64l) have no solution, so 
at these values of the cavity degree K all energy levels have zero width in the bulk of the 
sample. At a; = 1 the system of equations (I65II66I) can be solved explicitly: 

K*{g) = 2ge^l^^''^ (67) 
z*{g) = 2eg (68) 
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Notice that for small g the value K*{g) is exponentially smaller than the value Kc{g), given 
by f|T8|) . at which superconductivity appears. 

Next, we consider the region of small u. To solve the system ( |65]|66|) in this limit, we 
first notice that at a; = the solution of the equations gives the critical point discussed in 
Sect Jill PI -^'3(0) = Kc{g) = ge^^^''^'^ and ^(0) = eg. At small w < 1 and g <^ I the exponent 
z should remain small: z <^ 1 which allows us to neglect terms linear in u in fl65|66p : 

e(^/e9)-i [1 + (^/2)-] = ± (^] ' (69) 

eg \ K J 

~ ln(u;/2) = l-A + ^ln§ (70) 



1 + {uj/2Y ' ' eg K 

The second equation flTOj) shows that (1 — z/eg) scales as cu^, so deviations of z from its 
critical value can be neglected in the first equation where terms linear in (1 — z/eg) cancel. 
This gives 



^d{9, K) 



Kc{g) 



eg 



- 1 



(71) 



K 

This reduces to a power law behavior of u!d{g, K) when K is close to Kc and Ud ^ 1: 

/ \ 1/(69) 

c.,(^7,i^) = 2(e^)V(^^)(^l--^J (72) 

These results are illustrated by Fig. 8 which shows u)d{g = 0.129, i^) as function of K. 

Note that the previous analysis shows that the region of small w — 2^ gives a negligible 
contribution to /^j. This gives an a posteriori justification to the fact that we have neglected 
the non-perturbative modifications of the equation fl60|l in our study of the insulating phase. 



B. Scaling of the level width close to the transition. 

The level widths strongly depend on the frequency for uj > uJd{g)- The form of this 
dependence is important for the study of low temperature properties discussed below. To 
find T{uj) we need to rederive the equation ( |60l) for the case of non- negligible r(ci;). The 
computation using the Keldysh formalism, which we shall explain in Section IIVCI below, 
gives at zero temperature: 

T. = i2g/Krj: (73) 
k(i) - + ^ k 
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FIG. 8: Critical frequency, w^j that separates zero-width levels from the levels with non-zero width, 
plotted as a function of K for g = 0.129. This value of g was chosen because it corresponds to 
the experimentally relevant transition temperature at large K: Tbcs = 0.001 (see section rvT|) . At 
K < K* all states have zero width. In the intermediate regime K* < K < Kc states to the left of 
Ud{K) line have zero width (infinite decay time). The insert shows a blow-up region around Kc- 



As one expects, the decay rate of the state k provides the cutoff of the divergence at w — )■ 2^k- 
This equation is similar to the equation ([7]) for the fields appearing in the superconducting 
phase. We are interested in the scaling of the typical level width for u > Ud{g, K) but close 
to the transition. 

We will analyze this scaling with the same method as in Sect. IIII Ft in the RSB region. 
The crucial ingredient is the distribution of widths iy(F). We shall assume that it has a 
power law form with an upper cutoff : 

J- n 



W{T) 



where Fq < F < F^ax 



The non-linear mapping (j73|) gives a self-consistent equation for this distribution. Performing 
the same steps as those used in the derivation of Eq. (Hi|) . we get: 



9 



2b 



1 



where the function 'j{x,y), defined by: 



dT W(T)T^ 



K (2g_ 
K 



2b 



ciFiy(F)FS(26,^) (74) 



dt 



(75) 
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reduces to 7(x) when ?/ — )■ 0. The mapping (173|) imphes that values F much larger than g/K 
are very rare. Indeed, such a value of F at one step induces a much smaller value at the next 
step. We thus assume Fmax ^ g/ K ^ u^- Evaluating both sides of Eq. (17i|l we find 



2b 

9 



1 



^ K f2g\"' 7(2&) . ^ 



Close to the critical point K = Kd{g), the exponent 2b ^ m = 1 — eg. Expanding 
{g/gd{u!, K))"^ — 1 at small {u — ujd)/ijJd we find using the Eq.( l72|) that 



^gd{uj,K)J \ KcJ uJdig,K) 

Finally, we obtain for the typical level width: 

r^y^iu) ^ Fo(u;) ^ e^'/^^ exp f ^^^l (78) 

As we shall discuss below, this fast energy dependence of the level width has important con- 
sequences for the low-temperature transport properties. As before, the numerical coefficient 
C ~ 1 cannot be determined analytically because we do not know the precise value of Fmax- 
We emphasize that, according to Eq.( ff9l) . the energy scale ui in the exponent is much larger 
than the threshold energy Ud- This inequality uji ^ ujd is valid in the whole range of validity 
of Eq. ( ff9l) . as long as < g/K. Qualitatively, it leads to a very sharp growth of the 
typical level width Tq{ijj) right above the threshold. 

The estimate of the level width close to the quantum critical point, Kc{g), requires a 
special treatment because the expansion ( iTTll and the final result (|78ll are valid only for 
frequencies u close to the threshold Ud so that u — Ud ^ oud- Thus, they are not apphcable 
at the critical point K = Kc{g) where both Ud and ui vanish. An approximate expression 
for Tq{uj) in this regime should be found differently. First, using the function u)d{g, K) given 
by Eq. fl72|) . we determine the inverse function gd{uj, K) for low u: 

gd{uj,K) fuya- 1 



gc V 2 / 1-egc 

Next, we substitute this expression into Eq. (!76l) and find 

where T{gc) is some function of gc only. At egc ^ 1 we find T{gc) = C{e^gc)~^. 



To{uj) = ^exp 



-rigc) - 



(80) 
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C. The effect of a non-zero temperature on the level width. 



In this section we derive the cavity equations for the level width within the Keldysh 
formalism. This allows to study the effect of a non-zero but low temperature. A low 
temperature affects the relaxation rate in several ways. First of all, it changes the occupation 
numbers of the excited states and the ground states; this affects the perturbative equation 
fl60|) and thus shifts the position of the Ud{g) line. This effect is however small at T ^ 1. 
A more important effect is that, in the intermediate phase K*{g) < K < Kc{g), a non-zero 
temperature induces a small number of mobile excitations with frequencies above Udig, K). 
These excitations provide a mechanism for a small but non-zero level broadening for the 
levels even at very low frequencies, u < uJd{g)- 

The effect of a small number of mobile excitations can be estimated qualitatively by the 
following arguments. The excitation with energy E moves with a typical rate r*^^(co') ~ 
exp{—ui/{E — Ud)) from site to site. Thus, with exponential accuracy a typical site sees 
a mobile excitation with energy E passing by with rate exp{—u)i/{E — uod) — E/T). The 
dominant contribution to the relaxation of this site comes from the energies E = Ud{g) + 
y/ooiT. It results in the temperature dependence T ~ exp(— 2a/cji/T — cod/T) that shows a 
crossover between a behavior exp(— C/a/T) and an Arrhenius behavior exp(— C/T) as one 
goes away from the critical point. 

We shall now use the Keldysh formalism to support these qualitative reasoning. 

1. Cavity equations for the relaxation rate: derivation with the Keldysh formalism. 

As a first step, to be used later in the cavity approach, we begin with a study of a reduced 
two-spin system where a fluctuating field h{t) is coupled to spin cti. The Hamiltonian 
describing this system is 



We need to find the effective relaxation rate of spin 0. For this we employ the Jordan- Wigner 
transformation, using Fermion creation operators Cq and c\: 



H 



.+ 
1 



(81) 




4(1 - 2clci), (T^ = 24co - 1 
c|, al = 2c\ci - 1 



(82) 
(83) 



a 
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The Keldysh action obtained after averaging over the environment variables is: 

S = ^ {coaidtCoa + CiaidtCia - aH[c])dt / Cap{t -t')Cio,{t)cif}{t')dtdt[M) 

H[Ca\ = -(2^oCOaCo„ + 2^iCic,Ci„ + {2g / K) {Ci^aCla + ClaCo„)) (85) 

Here subscripts = ± are Keldysh indices corresponding to fermions moving forth or 
back in time. It is convenient to rotate all fermion vectors in the Keldysh space according to 
the Larkin-Ovchinnikov prescription (for details of Keldysh formalism see the recent detailed 
review^). For each i G {0, 1} we define 

C^l{t) = ^(Q+(t) + Q_(t)) , Q2(t) = ^(c^+(t) - Q-(t)) (86) 

Qi(t) = 75(c.+ (t) - c,4t)) , c,2(t) = ^(c^+(^) + c^-{t)) (87) 
After this rotation, the action acquires the form 

^ = J2f i'^OatdtCoa + C^aldtCla " H[c])dt - J2 f ^^'^^ " t')Clait)Cuit')dtdt' (88) 

a a,b 

where a,b & (1, 2). The Fourier-transform of the matrix Cab{t — t') is a standard triangular 
matrix in Keldysh space: 

The functions C^{u), C^{ijj) and C^{co) are respectively the retarded, advanced and Keldysh 
component of the matrix C{uj). If we denote by ri(ci;) the spectral density of the external 
noise h{t) acting on spin 1, they are given by 

C^{uj) = -iTi{u) C^{uj) = iViiu) C^{u) = -2iri(w)tanh^ (90) 

Our goal is to determine the spectral density of noise acting on spin 0. The bare retarded 
fermionic Green function at site 1 (in the absence of the coupling to site 0) is 

G«(o)(a;) = (a;-26+^ri(a;))-^ (91) 

The action (l88|) is quadratic; diagonalizing it we get the imaginary part of the retarded 
fermionic Green function at site 0: 
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Here ei,2 = + ^2) ± y ~ '^2)^ + {'^g/KY are the eigenvalues of the two-spin Hamil- 
tonian (ISTl) in absence of the fluctuating field. In the following we shall neglect the effects 
of level repulsion which translates into the difference between actual energies ei, €2 and their 
unperturbed values ^1,^2- This is the same approximation that we employed in Sect. IIIII 
to obtain analytical results. As we have shown above, this approximation remains accurate 
even for very modest values of K. Using this approximation we can rewrite ( l92l) as 



as one expects from perturbation theory arguments. 

At zero temperature these results for the fermion Green function are translated directly 
into spin correlators. At non-zero temperature this conversion is less trivial because the spin 
correlators acquire an additional decay compared to fermions. Physically, this is due to the 
fact that the spin Hamiltonian is non-linear, so the thermal excitation of one spin might lead 
to the relaxation of another spin. For this process to happen the thermal excitation should 
be mobile, meaning that it should have an energy larger than Ud- This makes such processes 
rare at low temperatures. Formally, the retarded transverse spin Green function at site 0, 
denoted as Dq{u), differs from the retarded fermion Green function, Gq{u), because the 
former contains an additional factor Di =< (jfcr^ > due the non-linearity of the Jordan- 
Wigner transformation. More precisely, in Keldysh technique these two Green functions are 
related by 



The dynamics of the z-spin components is purely dissipative in the leading order in J 
which we consider here. In this approximation the retarded Green function of the z-spin 
components is zero; this allow us to neglect the second term in the general expression f lM|) 



for 0^(1). The first term has two parts: f-Df" = Tn? + |-Df where m =< >= tanh(^i/r) 
and f-Df^ is the irreducible part of the symmetrized spin-spin correlator at site 1: 




which implies that the imaginary part of the self-energy is 




(93) 




(94) 
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Using Eqs.dHa |9T]) we get: 



-D^{u) = 2(1 - m') ^^'^^^'^ , (95) 

2 ' ^ ^ ^ ^2 + [211(2^1)]' 



In the following we shall be mostly interested in the relaxation at low energies u ^ Ud- 
Combining Eqs. (l94|l95|) and ( 19T|) . we get the equation for the imaginary part of the retarded 
transverse spin-spin correlator at site at non-zero temperature: 



Here rQ°''(u;) is the relaxation rate at T = 0; it coincides with [G^(a;)] ^ given by Eq. (!93|l . 
Below we neglect the level widths compared to level energies, which allows to find the 
imaginary part of the inverse transverse spin correlation function at site 0, i.e. the relaxation 
rate of the spin at site 0: 

^ ° ' \kJ I (1^ - 26)' + r?(i^) (i^-26)' + 4r?(26)/ 

The second term in this equation is due to the non-linearity discussed above; in this term 
we have neglected the contribution of Ti{uj) compared to ri(2^i). The reason for this is 
that this term is exponentially small at low T due to the factor 1 — ^ 1. Thus, this 
term makes a significant contribution only if the energy of spin 1 is large: > ua , so that 
ri(2^i) is much larger than Ti{u). 

The result fl97|) describes the relaxation in the system of two spins. We now consider 
a full cavity problem in which spin i has K neighbors, labeled k = 1, . . . ,K, and each of 
these neighbors feels a fluctuating external field with spectral density Tkioo)- We compute 
the relaxation rate of the spin at site 0. Adding the contributions from all neighbors we get 

r (u) - V [ ^'r^(^) + (i-m2)2r,(2G) 1 , . 

^\(c.-26)Vr?(c.) + (c.- 26)^4^(2^0 J ^ ^ 

In the zero temperature limit ruk — >■ 1 the second term in (|98|1 vanishes and this formula 
reduces to ( [73|) . 

At non-zero temperatures, one needs to take into account the second term in fl98|l . which 
becomes important for low frequencies u < Ud for which Tk{uj) would be zero at T = 0. In 
this regime the second term works as a source term to the recursion that would otherwise 
give zero. 
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2. Solution of the recursive equations with a source term and consequences for low temperature 
properties. 



At non-zero temperature the low energy modes which are discrete at T = acquire a 
finite hfetime. We shall estimate their broadening at very low temperature when the effect 
of non-zero T on modes with u > Ud can be neglected. Our starting point is equation f l98p 
which contains two terms of very different physical meaning. The first term describes the 
decay due to the indirect coupling to the external fields far away, it is the only term present 
at T = 0. The second term describes the relaxation caused by mobile thermal excitations of a 
neighboring spin with energy above the threshold Ud- At low temperatures such excitations 
occur exponentially rarely and thus are essentially limited to a narrow range of energies 
slightly above Ud- Because the density of these excitations is very low, we can ignore the 
situations in which a given site feels more than one thermal excitation nearby. It implies 
that the relaxation rate induced by these excitations for low-energy modes is much smaller 
than the rate of the mobile excitations. This allows us to replace the effect of the mobile 
excitations on low energy modes by an effective fluctuating field. As we shall see below the 
spectrum of this field is featureless, so in all respects the mobile excitations are similar to 
the external field. 

The approximation discussed above allows us to replace the relaxation rate rk{2C,k) in the 
second term in fl98l) by its typical value and to ignore the effect of a non-zero temperature 
on this term. Eq. (!98l) becomes very similar to Eq. (l48l) for the susceptibility: 



kii) - ^^k) + i k 



where Ttyp{x) is given by the zero-temperature result, see Eq. (ff8l) . The dominant contribu- 
tion to the external dissipation ?7(T, uo) in Eq.( !99l) is determined by the competition between 
exp(— s/T) and Vtyp{x) ~ e~^^^^ {x — Ud)]- Evaluating the integral over x in the 

saddle point approximation (valid under the conditions uo < Ud and T ^ wi ) we find that 
ri{T,u) is weakly w-dependent and given by 



T V T 



(100) 
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The estimate fllOOp is valid with exponential accuracy provided that T ^ wi. The weakness 
of the frequency dependence in ri{T,u) ~ vi^) proves our statement above according to 
which the noise produced by the thermal excitations is featureless and exponentially small. 

Using the formal analogy between Eq. fl99l) and the Eq. fj49|) for the susceptibility, we 
can now use the results obtained in Sec JIlIHl The key conclusion of this section that we 
need here is that the typical value of the susceptibility does not contain any divergence 
at the T = transition point. This implies that the typical sub-threshold level width is 
r*^^((X' < Ud) ~ vC^)- This confirms our conjecture above that the relaxation rate of the low 
energy modes is very low, which allowed us to replace the mobile excitations by an effective 
fluctuating field. To estimate the level width for the whole range of frequencies, we note 
that both ri(T) and ro(a;) (see Eq. flTSl) ) are exponentially fast functions, therefore 

^ max (?7(r), ro(a;)) . (101) 

The estimates ( j78l) and ( llOOp are not valid right at the critical point g = gc, where Ud and ui 
vanish and the level width is given by Eq. (l80!) . Repeating the calculations similar to those 
used to derive Eq. (l80l) we find at g <t^l, instead of (llOOp : 

(102) 

The equations (110011102^ give the typical relaxation rate of low energy modes. At the 
same time they give a typical noise level and thus a typical rate of transport processes in 
this model. For the superconductor- insulator transition the result ( llOOp implies that, away 
from the critical point, the resistivity follows the Arrhenius law at very low temperatures 
and a exp(l/\/T) law in the intermediate temperature regime. This behavior is exactly 
opposite to the one expected and observed in conventional Mott insulators where Arrhenius 
is followed by Mott behavior as temperature is decreased. Exactly at the critical point, the 
resistivity grows as a stretched exponential exp(T^") with a ~ 0.25 — 0.3 for a physically 
relevant g ~ 0.12. 

V. THE EFFECT OF FRUSTRATION. 



TjiT) oc exp 



C 
e^g 



Disordered superconducting films close to the SI transition can be driven to insulators 
by the application of a magnetic field. The properties of these materials in presence of a 
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magnetic field have been extensively studied experimentally and are rather unusual. So it is 
important to discuss the theoretical expectations. We will not undertake a full quantitative 
study here, but we will only study the leading effects of a magnetic field on the phase diagram 
of our model formulated on the Bethe lattice. In the framework of this model the effect of 
a magnetic field is described by the effective model where nearest-neighbor couplings Mij 
defined in Eq.([T]) acquire random phases: Mij = ■^je^"', with (a^) = /. The effects of these 
random phases on the two major lines of our phase diagram (see Figs. [2] and 8), the critical 
temperature line Tc{K) and the threshold energy line in the insulating phase uJd{K), are 
crucially different. Whereas all the equations for level widths Fj contain squares of absolute 
values of matrix elements |Mjjp only, and thus do not depend on /, the equations for the 
order parameter are affected by the random phases. 

Consider first the case of small phase fluctuations / ^ 1. In the limit of very large 
K > K^^^{g) the simple mean-field approximation should be valid, and the transition tem- 
perature is determined by the equation for the first moment of the P{B) distribution. Ran- 
dom phases enter this equation via a straightforward modification of the coupling strength, 
g — )> g{l — //2), which leads to the suppression of T^. 

\^IA11 = -L (103) 

Tc 2g ^ ' 

Because g is small, the suppression of the transition temperature (and therefore the decrease 
of the typical order parameter in the ordered phase) can be rather strong even at / ^ 1. 

The case of strong frustration is more complicated. Here we consider the limit of very 
large fields that generate completely random phases ai with uniform distribution over (0, 2tt). 
Instead of Eq. (j7]) we get 

S-^lE^^Tf'-h/Jv'^. (104) 

k=i y^k + ^k 

and we look for a solution for the P{B) distribution function which depends only on the 
absolute value \B\. 

In order to determine we use the linearized version of f ll04p which can be rewritten in 
terms of the Fourier transformed distribution Q{s^s*) = J dBdB* P{B^B*)e^^^^*^^*^\ We 
shall assume that this Fourier transform depends only on the absolute value Q{s, s*) = 
Q{\s\). It then satisfies: 

. g tanh "' ^ 
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(105) 



Note that these equations are formally identical to those obtained in the case without mag- 
netic field in f l29|) . However, the analytic properties expected in the present case are different: 
the random phases induce a symmetric P{B) distribution, instead of the distribution sup- 
ported on positive 5's when there is no magnetic field. We will look for a power-law solution 
of Eq. fllOSp at small s in the form Q{s) = 1 — 

Because the average order parameter is zero in the limit of large magnetic fields, the simple 
mean-field solution is obtained by assuming that P{B) is determined by its second moment 
and correspondingly x = 2. The self-consistent equation f llOSp then gives T^^^ = 1.705g^/K. 
This is the equivalent of @ in the zero magnetic field case. 

The equation for the Laplace transform (11051) is formally identical to the one obtained in 
zero field f l29|) . so the value of the exponent x that determines the behavior of the Fourier 
transform at small s is determined by the same equations fl23]IM|) derived in Sect. UTTDI 
The important difference between the cases of zero and large magnetic field is due to the 
fact that in the latter case the simple mean field solution corresponds to x = 2. As a result, 
in contrast to zero field case, the simple mean field solution is not valid even in the K ^ oo 
limit. In other terms, RSB always occurs for the fully-random problem defined by Eq. ( 1104^ . 

We begin by solving these equations in the large K limit, where the exponent x approaches 
1 + e, with e ^ 1 as long as g <^1. Assuming that eln ^ ^ 1, we can extend the integrals 
over ^/T in Eqs. fl23|24l) up to oo. We estimate the resulting integrals: 



Similarly to the zero magnetic field case, in the RSB phase the naive mean field prediction 
T^^ = 1.705g^/K is exponentially larger than the correct result f ll07p for small g. 

We now prove that for any K the transition occurs in a RSB phase. To find the tem- 
perature, TfisB, corresponding to replica symmetry breaking we consider the Eq. (1241) and 
assume that x = 2, using a procedure similar to the determination of the RSB point, Eq. fl2T|l . 




and obtain 



e 



(106) 



d-l/eg 

K 
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discussed in Sec JIIIDl for the unfrustrated case. However, in contrast to the zero-field case, 
the corresponding temperature TfjSB ~ g/K, is always above the simple mean field value 
T*^-^ at g 1. Thus, it is necessary to solve both equations fl23f24p together to determine 
the value of the exponent x < 2 and transition temperature. 

The applicability of the solution (I106II107P is limited to the regime of very large K, when 
T is so small that the corrections of order T*^^ due to the finite upper limit in the integral 
lo^'^ ( ^ ) ^^^^ cire negligible. When K decreases these corrections become significant and 
the exponent x starts to decrease; it eventually approaches unity at the value K = K^^^ 
determined in Eq. fl2T]) . At the same time, Tc{K) deviates from the simple law fll06p and 
attains at = K^^^ its maximum value, equal to the mean-field transition temperature 
Tco of the unfrustrated model. At still smaller the solution for T^. is identical to the one 
discussed in Sec JIIIDl for the unfrustrated model. This behavior is summarized in Figj9l 
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FIG. 9: Blue line: dependence of the critical temperature Tc{K) on K of the fully frustrated system 
for g = 0.1. Red line: result of the simple mean field approximation for the problem without 
frustration. The maximum of the blue line corresponds to the case where the order parameter 
distribution decays with an exponent x = 1 and K = K^^^, where K^^^ is determined for the 
unfrustrated model by Eq. (|2ip . At K < K^^^ the critical temperatures of the regular and the 
frustrated model coincide. 

Although the transition line Tc{K) stays the same (in the RSB phase) for the unfrustrated 
model and for the strongly frustrated one, the amplitudes of the order parameter in the 
ordered phase differ considerably. This is illustrated by the numerical solution of the random- 
phase mapping equations (11040 shown in Figs. [10] and [11] for zero temperature and near to 
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the quantum critical point g = Qc- 
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FIG. 10: Distribution function of the local order parameter for the strongly frustrated model with 
K = 4 at three different values of the coupling constant. The number near each curve indicates 
the corresponding value of {g/gc — !)• The plateau in P{B)B^^^ demonstrates the existence of 
the scaling regime ()3ip . 
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FIG. 11: Typical amplitude of the order parameter as function of the proximity to the quantum 
critical point. 

The results shown in Fig. [11] are very much like those seen in Fig. [3] for the unfrustrated 
model, but the typical ampUtude B^yp is suppressed by a factor ~ 100 in the frustrated case 
at the same g/gc value. 

To summarize this section, we have demonstrated that frustration suppresses strongly the 
transition temperature when one is in the replica symmetric phase at sufficiently large 
K but it has no effect on in the RSB phase near the quantum phase transition. This 
can be interpreted as a consequence of the quasi-one-dimensional nature of Bethe-lattice 
clusters which contribute to the formation of the coherent state in the RSB phase. While 
Tc is unchanged in the RSB region, the amplitude of the order parameter (at T ^ T^) is 
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strongly suppressed due to frustration, as demonstrated in Fig. [TTl More work is needed to 
decide if these results, obtained on the Bethe-lattice problem, can be applied to the finite- 
dimensional problem where closed loops are present (see discussion in section IVl Ch We 
expect, however, that the results will remain qualitatively similar due to the dominance of 
a small number of paths that makes the presence of small loops largely irrelevant. 

VI. CONSEQUENCES FOR EXPERIMENTS. 

A. Distribution of coherence-peak heights in STM and Andreev point contact 
tunneling experiments. 

One of the main results of this work is anomalous broadening of the distribution of the 
local values of the order parameters in the vicinity of the SI transition. This conclusion can 
be tested by STM measurements. 

The same STM experiments can also confirm that the transition happens due to Cooper 
pair localization and not by unbinding of Cooper pairs. In a usual setup of these experi- 
ments one measures the current- voltage characteristic of a highly resistive tunneling contact 
between the material and a needle of the STM machine. The measured differential con- 
ductance, dl/dV, is proportional to the single-electron density of states, p{E), at energy 
E = eV. Thus, these measurement directly probe the superconducting gap. The exper- 
iments performed on superconductor-insulator transition in TiN— and InO^^ demon- 
strate that in these materials the superconductor-insulator transition happens without pair 
destruction (the single-electron gap remains large). The theoretical justification of this was 
given in the work^S and is summarized above in section III A[ Thus, one expects that this 
transition is driven by the competition between disorder and Cooper-pair tunneling and is 
described by the models studied in this paper. As we discussed in section Hi Al the present 
paper ignores subtle effects of the correlations between matrix elements, in this approxi- 
mation the disorder changes only the number of neighbors, K. In order to compare with 
experimental data we need to choose the interaction constant so that to reproduce the cor- 
rect value of Tc/Ep ~ 10~^ away from the transition in these materials. This gives the 
interaction g = 0.129 used in many numerical plots in the previous sections. 

In superconductors the density of states above the gap displays the coherence peak which 
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is due to the superconducting order parameter. The height's statistics of these peaks can be 
used to measure the variations of the local order parameter appearing in the vicinity of the 
superconductor-insulator transition. In particular, the broadening of the P{B) distribution 
translates into the broad distribution of peak heights. This theoretical prediction was indeed 
confirmed in very recent experiments^^ where superconducting samples with different values 
of Tc were studied and compared. The data shows that while the superconducting gap 
remains largely intact in lower Tc samples and experience modest fluctuations from point to 
point, the distribution of peak heights changes dramatically as the insulator is approached. 
This observation is in full agreement with the theory developed in this work, and confirms our 
expectation that the solution of the model on Bethe lattice provides a good approximation 
for the properties of realistic physical systems. 

Another possible set of data is given by Andreev point-contact tunneling measurements. 
Similar to conventional tunneling, these experiments give current-voltage characteristics of 
the point contact. However, in contrast to conventional contacts, the contacts used in these 
experiments are characterized by a small resistance, R ~ Ikfl. This allows a coherent 
tunneling of the Cooper pair between the normal needle and a superconducting material. 

In a conventional BCS superconductor the differential conductivity observed in Andreev 
point-contact experiments shows a suppression below the single-particle gap and a peak at 
the gap edge, similar to the conventional tunneling. The reason for this is that in conventional 
superconductors the gap, 6Epair, for collective pair excitations (pseudospins) is exactly twice 
as large as the gap for single particle excitations: SEpair = 2Abcs- Tunneling of Cooper 
pair brings energy 2eV, so the gap edge observed in differential conductivity for the single 
particle tunneling coincide with the one observed in Cooper pair tunneling. 

In contrast to conventional superconductors, strongly disordered superconductors in the 
vicinity of the SI transition are expected to show much smaller gap for collective pair exci- 
tations, SEpair, than the one characterizing single-particle tunneling, Agp. In this case, the 
observed differential conductivity below single-particle gap is due to the coherent Cooper pair 
tunneling which might occur at potential differences eV < Agp. In the lowest non-zero order 
in tunneling one gets that at such bias the differential conductivity is proportional to the 
local susceptibility, lmxr{^ = eV^), of the equivalent spin model evaluated at the position of 
the needle, r. This expectation is correct provided that the collective modes contributing to 
the local susceptibility are delocalized at this frequency. The tunneling into localized modes 



49 



cannot lead to experimentally observable conductivity. Because the energy separating the 
localized and delocalized modes is the same at different positions in the sample, we expect 
that the observed differential conductivity is characterized by a spatially uniform threshold 
or a peak. This expectation is in agreement with the preliminary results^. 

B. Low-temperature resistivity in the insulating state 

Because STM measurements require that the resistance of the tunneling contact is much 
larger than the resistance of the sample, they become difficult, if not impossible, in the 
insulating phase. So, in this state one should compare the theoretical predictions with 
the transport data which are always more difficult to interpret. The most striking feature 
of the resistivity on the insulating side of the transition in InO and TiN is the activated 
behavior of the resistance R ~ exp(T4/T), with an activation energy that decreases in the 
vicinity of the transition.—*^ This decrease of the activation energy, T^, is in contrast with 
a large and non-critical single particle gap observed in STM— li^ii^ on the superconducting 
side of the transition. It indicates that the transport in this regime is due to incoherent 
Cooper pair hopping. More direct evidence of the pair transport is provided by the data on 
magnetoresistance oscillations in perforated films which displays periodic oscillations with 
"superconducting" flux period $o = hc/2e in the insulating state.-^ 

As we have shown in Section IIV At the low lying excitations of the spin model have 
zero width. This translates into the localization of all excitations with energies E < Ud- 
The transport is made possible by modes with energies E > Ud which become thermally 
excited at non-zero temperatures. Thus, we expect that the resistivity in this state has an 
approximately activated temperature dependence at lowest temperatures with an activation 
energy = ^d which decreases when one approaches the transition. The subdominant 
term in the resistance is proportional to exp ^Jui/T which might become important at 
intermediate temperatures because uji ^ ujd- This result is in agreement with the data. It 
should be contrasted with the naive expectation according to which the transport of Cooper 
pairs could be characterized by Mott (or Efros-Shklovskii) behavior, giving a behavior R ~ 
exp(TA//T)'* (a = 0.25 — 0.5), in analogy with single electron transport in conventional 
insulators. In these insulators the activation behavior is realized at intermediate temperature 
range and Mott behavior at lowest temperatures. 
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The sharp boundary, w^, separating the locahzed and delocahzed modes could be more 
directly probed by microwave experiments. We expect that close to this boundary the 
lifetime of the excitations becomes very long. Thus, even a weak microwave radiation could 
excite a large number of these excitations. This would dramatically increase the conductivity 
of the sample without affecting significantly its temperature. Furthermore, applying the 
microwave power in pulses with varying time intervals between the pulses, one could obtain 
information on the lifetime of these excitations near the many-body mobility edge. 

C. The effect of frustration induced by a magnetic field or by unconventional order 
parameter symmetry. 

As discussed in Section IVj the highly inhomogeneous spatial structure appearing in the 
vicinity of the SI transition strongly suppresses the effects of frustration and thus of a 
magnetic field. In the approximation used in this work the magnetic field has no effect 
on the transition temperature in this phase. The account of higher order effects m. 1/K 
would lead to a small suppression of the critical temperature in this regime. This conclusion 
is in qualitative agreement with the data: it has been found (see e.g.—"— and references 
therein) that the properties of the materials depend very sensitively on the disorder level 
in the vicinity of the SI transition in contrast to their smooth dependence on a magnetic 
field. In spite of its weak effect on the transition temperature, we have observed that the 
magnetic field should affect strongly the order parameter. This prediction can be tested 
by STM and Andreev point-contact measurements which provide direct information on the 
local superconducting order parameter and its distribution. 

The very weak effect of a magnetic field on the transition temperature implies that su- 
perconductivity might survive in strongly disordered superconductors with unconventional 
pairing symmetries such as d-wave symmetry, provided that the Cooper pairs in these super- 
conductors are formed on very short scales which are weakly affected by disorder. This might 
be the case in high T^. superconductors in the pseudogapped regime in which there are rea- 
sons to expect^*^ that electrons are paired in the corners of the Brillouin zone. In this case 
the Cooper pairs have very small size and are not affected by disorder. The global supercon- 
ductivity is due to the coupling between these localized Cooper pairs. These couplings have 
random signs due to the d-wave symmetry of the order parameter which makes this situation 
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similar to the strongly disordered superconductor in a magnetic field. As discussed above 
the randomness of the couplings has weak effect on the transition temperature. This might 
explain the apparent similarity between underdoped high superconductors and strongly 
disordered InO and TiN films, noted by many researchers previously, see for example^i^ 

D. Change of level statistics: suggestion for numerical work. 

The many-body mobility edge line Ud{K,g) found in Section H V Al divides the parameter 
space of the Hamiltonian into regions with qualitatively different spectral statistics. At 
u < Ud the spectrum is point-like and the many-body wavefunctions are localized, thus 
we expect the level statistics to be Poissonian in this region. On the other hand, above 
the mobility edge, at w > Ud, a nonzero line width is generated, indicating the extended 
nature of wavefunctions. In this regime we expect a Wigner-Dyson level statistics, with level 
repulsion. This change of level statistics can be studied numerically via exact diagonalization 
of the Hamiltonian ([1]) . 

A similar study was implemented recently in the work^^ for a one-dimensional model 
of interacting fermions. This paper has found a qualitative evolution between Poissonian 
and Wigner-Dyson statistics, but large finite-size effects make any quantitative conclusions 
difficult in the absence of an appropriate analytic theory. We expect that the Hamiltonian 
([1]) with moderate branching number K > 1 could be useful for such study, because our 
theoretical predictions for the position of the mobility edge can be used for comparison with 
the numerical data. 

As a next step, it would be important to extend such a numerical study to the same 
type of quantum spin models on a usual real-space lattice. The comparison of the results 
for spectral statistics obtained on conventional Euclidean lattices and on the Bethe lattice 
could be very useful to check the applicability of the approach developed in this paper to 
physical systems. 

VII. CONCLUSIONS 

In conclusion, we have outlined a solution of the strongly disordered spin model on the 
Bethe lattice that can be mapped onto the disorder-driven superconductor-insulator tran- 
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sition or onto the disordered ferromagnets. We found a series of two zero temperature 
transitions between an ordered state, a state with a slow relaxation and a state with no 
relaxation. Our solution also shows that the low temperature phase in this model is always 
very strongly non-uniform with both the order parameter formation and the spin relaxation 
controlled by rare paths that contain a very small fraction of spins. When applied to the 
superconductor-insulator transition our results imply the existence of both weak and strong 
insulators. In the former the relaxation rate varies faster, as exp{y^7Ji/T) , in an intermediate 
temperature range, but it crosses over to Arrhenius exp{u)d/T) at very low temperatures. 
Exactly at the quantum critical point we expect a relaxation rate that varies as exp(l/T)" 
with a non-universal a that depends on the interaction constant. For physically relevant 
values of the interaction constant we expect a ~ 0.25 — 0.3. 

In the strong insulator the relaxation is completely suppressed. Of course, the physical 
effects neglected in our model (like electron-phonon coupling) would lead to some slow 
relaxation even in the latter phase. 

The predictions of our theory can be tested by conventional STM experiments and An- 
dreev point-contact tunneling experiments on superconducting films of InO and TiN as 
discussed in detail in Section IVll 

In the strongly insulating phase the local quantum levels remain discrete, this implies 
that such systems would remain coherent for a very long time (limited by the interaction 
with phonons and other environment neglected in our model). When applied to disordered 
magnets this result implies that a stronger disorder might make the spin system less noisy 
and more coherent. This might be important for the efforts to suppress the flux noise in 
superconducting circuits where it is believed^i^ to be due to the disordered spins at the 
surface of the superconductors. 

Our theory is directly applicable to the models on Bethe lattices which do not have small 
loops. Qualitatively it looks very likely that the dominance of a small number of paths, 
which is a central result of this work, implies that small loops are irrelevant and therefore 
the Bethe-lattice approximation should be very good. This conjecture might be checked 
numerically (see Section [VI Pp . 

We now briefly discuss possible extensions of the present theory and list few open prob- 
lems: 

• Extension of the developed technique to the spin lattices in Euclidean space. It is be- 
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lieved that the exponential behavior of physical properties in the localization problem 
on a Bethe lattice disappears in any finite dimension^. However, the arguments 
are not directly applicable to the non-linear problem of phase transition that we con- 
sider here in which the phase-space dimension increases exponentially with the number 
of sites. Qualitatively it seems likely that the dominance of a small number of paths 
is an indication that the main results of the Bethe-lattice approximation will remain 
correct in finite dimensions. 

The extension to finite-dimensional systems would be greatly simplified if one could 
find the appropriate order parameters for the phase transitions found on Bethe lat- 
tice. In particular, these order parameters should describe not only the appearance 
of the superconducting order but also its homogeneity and its time-reversal symmetry 
breaking which appears at the second transition. 

Direct calculation of transport quantities such as thermal conductivity or electrical re- 
sistivity near the quantum transition studied in this work. The present paper computes 
the spin relaxation time which for the SI problem translates into the time required for 
the Cooper pair to move away from a given site. With exponential accuracy it should 
be the same as the resistivity or thermal conductivity in the system. In order to com- 
pute the properties with better accuracy one would need to develop the formalism to 
compute transport properties directly in this problem. 

Direct computation of the dynamic spin susceptibility Xr{^)- The behavior of this 
quantity determines the pair conductivity measured in STM experiments and Andreev 
contact spectroscopy. 

Finite temperature effects in the strong insulator. In the approximation used in this 
paper even a high temperature has no effect on the strong insulator. However, one 
might worry that our approximation might miss exponentially small effects such as 
rare regions characterized by larger density of states that would make this part of the 
system less insulating. 
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IX. APPENDIX A: QUANTUM CAVITY METHOD: SUZUKI-TROTTER FOR- 
MULATION AND FURTHER APPROXIMATIONS 

If one is interested only in the thermodynamic properties of the problem described by 
f|T|2|) on the Bethe lattice (a random graph of connectivity Z = K + 1), one can write^^ 
an exact (at the replica symmetric level) cavity recursion in terms of Suzuki- Trotter paths. 
Unfortunately this exact recursion can be studied only numerically and is rather heavy to 
handle in the case of disordered systems. Our approach uses a simplified version which 
amounts to using a special measure on the Suzuki- Trotter paths. In this appendix we briefly 
summarize the exact cavity recursion and we make explicit the extra approximations of our 
method. Let us study for instance the transverse- field Ising Hamiltonian ([2]). 

Hi = -J2 ^^^^ - ^ E ^^^^ (108) 

The partition function Z = Tr e~^^'can be expressed with the Suzuki- Trotter repre- 
sentation. Using Nt imaginary time steps, one introduces the time trajectory of each spin, 
crj(t) G ±1 , where t = 1, . . . ,Nt. Then 

Z = \im J2 (109) 

Nt^00{(7,{t)} 

where: 

= ^^^(^)^^w - EEr^^^(^)^^(^+ 1) (110) 

* t (ij) t i 

and Ti = i log cosh (^^) 

The exact RS method introduced in^"^ is the following: take a branch of the Bethe lattice 
rooted on spin i, and define the probability distribution of the time-trajectory of this spin as 
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ipi [(Ti{t)]. Let us define analogously, for each of the K spins j which are the first neighbors 
of i on the rooted tree, their time trajectory as tpj [aj{t)]. Then one can write formally the 
mapping that generates, from {ipj [cTj(t)]} the new ipi [<Ji(t)]. This mapping is then typically 
studied numerically by a population dynamics method where one represents each of the 
ipj, V^j by a sample of spin trajectories. In paper'^^ it is shown how to take the Nt = oo limit 
by using a continuous time formulation where one memorizes only the times at which the spin 
trajectories jump. The whole procedure can be seen as a two-step approach: transforming 
the problem to one of classical spin trajectories, and then using the classical cavity method^^. 
In particular, the cavity mapping can be shown to derive from a variational principle based 
on a Bethe free energy. 

The quantum cavity approach that we have introduced in Sect. Illl Bl can be understood 
as an approximate way to study the exact RS Suzuki- Trotter based cavity method. Basically, 
the hypothesis amounting to using a local Hamiltonian of spin k in the form C.k^'l + B^a^ 
amounts, in the Suzuki- Trotter formalism, to using a local distribution ipk [crfc(t)] which takes 
the special form 



This form can be used with different methods. In a variational approach one injects this form 
into the Bethe free energy. This gives a free energy which depends on all the parameters Bk, 
on which one should optimize. In practice this method is a bit heavy numerically because 
the computation of this Bethe free energy as function of the Bk involves, for each spin i, a 
trace over aj and all its neighbors 0"^, which is thus a sum over 2^"*"^ terms. In the approach 
that we have described in Sect. IIIIBI we use a slightly different method: we suppose that all 
the cavity neighbors of a given spin i have distribution tp^ Wk{t)]- From these distributions 
we compute (af ) and deduce from it the value Bi which reproduces this average. This 
defines a mapping from {B^} to Bi. This type of approach can be tested by comparing 
it to the full numerical sampling of Suzuki- Trotter trajectories in the case of a uniform 
transverse field where = ^. This study, done ib^, confirms that this approach is able to 
reproduce the phase diagram accurately (the variational approach is a bit more precise, but 
both approaches give the zero-temperature value of Qc within a few per-cent). 

Our approach uses one more step of approximation: instead of the quantum cavity map- 
ping described above, we have mostly used the explicit mapping ([7]) which has the advantage 



^k [cTfc it)] = C exp 




^ J2 + /^r. E akit)akit + 1) 



t t 



(111) 
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that one does not need to compute (erf) in order to find the value of Bi. The vahdity of this 
further step of approximation has been studied in Sect JIII Gl 
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